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Motivated by recent experiments on quantum mechanical charge pumping in a Cooper pair sluice, 
we present a measurement scheme for observing shifts of transition frequencies in two-level quantum 
systems induced by broadband environmental fluctuations. In contrast to quantum optical and 
related set-ups based on cavities, the impact of a thermal phase reservoir is considered. A thorough 
analysis of Lamb and Stark shifts within weak-coupling master equations is complemented by non- 
perturbative results for the model of an exactly solvable harmonic system. The experimental protocol 
to measure the Lamb shift in experimentally feasible superconducting circuits is analysed in detail 
and supported by numerical simulations. 
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I. INTRODUCTION 

Realistic quantum systems are never completely iso- 
lated. Even a single atom in zero-temperature vacuum is 
influenced by the zero-point fluctuations of the electro- 
magnetic field which in turn induces a shift of its tran- 
sition frequencies known as the Lamb shift [l]. Cavity 
quantum electrodynamics (QED) provides a particularly 
convenient set-up to observe this shift since the restricted 
geometries of the cavities allow the atoms to interact only 
with the fluctuations of single harmonic fields [2 . Ac- 
cordingly, the exploration of fundamental quantum phe- 
nomena has reached an unprecedented level of accuracy 
and has extended our basic understanding of quantum 
mechanics [3l. Very recently, such a measurement has 
been performed in a solid state analogue, where a two- 
level quantum system implemented in the form of a su- 
perconducting Cooper pair box (CPB) is embedded in 
a superconducting waveguide resonator [4]. The advan- 
tage of this device is its tunability which allows to access 
ranges in the parameter space unavailable in quantum op- 
tical set-ups. The theoretical framework here is provided 
by the Jaynes-Cummings model and its recent extensions 
(see, e.g., piGl). 

In contrast to single-frequency environments, typical 
reservoirs for mesoscopic solid state devices are charac- 
terized by broadband spectral distributions in thermal 
equilibrium with no convenient access to manipulate the 
populations of individual modes. Energy exchange be- 
tween the system and the reservoir as well as other fluc- 
tuations lead to dephasing and relaxation, thus limit- 
ing time scales over which coherence and entanglement 
are preserved [7H9] . An understanding of these processes 
calls for the well-known theory of open quantum systems^ 
for which several approaches have been developed in the 
past. Among them, there are formally exact ones such as 
the path integral representation [9l [10] , and perturbative 
ones such as Redfield and master equations [7^. Within 



these weak-coupling formulations even explicit expres- 
sions for the reservoir-induced frequency shifts can be 
derived, while associated experimental observations are 
still missing. To fill this gap, in the sequel we discuss and 
analyse a theoretical proposal to retrieve the Lamb shift 
and its finite-temperature version (termed also ac-Stark 
shift according to the shift induced by a finite-amplitude 
drive) for a superconducting two-level system embedded 
in an ohmic environment. 

The two-level quantum system studied here is based 
on the Cooper pair sluice [11 , i.e., a superconducting 
island separated by two superconducting quantum inter- 
ference devices (SQUIDs) acting as tunable Josephson 
junctions. This device is embedded in a superconducting 
loop together with an additional large Josephson junc- 
tion employed as a threshold current detector, similar 
to the Quantronium device [12 . By adiabatically tun- 
ing the control fields of the sluice (left and right effective 
Josephson couplings and gate voltage of the island) along 
a closed path in the parameter space, the accumulated 
Berry phase [13] has been studied experimentally by ob- 
serving the transferred Cooper pair current as a function 
of the total phase across the sluice [14J. Furthermore, 
the impact of the environmental degrees of freedom have 
been investigated in a set-up where an additional resistor 
is coupled capacitively to the sluice island [15UT8]. The 
peculiarity of this architecture is that the temperature 
of the system is determined by that of the resistor. The 
same device is analysed here for a spectroscopic measure- 
ment; no time-dependent pumping signal is applied. This 
allows us to determine the Lamb and the Stark shifts for 
the circuit. 

This paper is organized as follows: In the following 
section, we briefly recall the weak-coupling approach for 
open quantum systems to introduce the notation and 
to give explicit expressions for the reservoir-induced fre- 
quency shifts. To validate our results and to gain insight 
in the regime of stronger system-bath interaction even 
without the commonly used rotating wave approximation 
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(RWA) , the exactly solvable case of harmonic degrees of 



freedom is discussed in Sec. |III[ Section [IV| presents the 
proposed experimental measurement scheme with numer- 
ical simulations specified in Sec.|V[ Conclusions are given 
in Sec. [Vll 



II. 



WEAK-COUPLING LIMIT AND LAMB 
SHIFT 

A. General theory 



The standard microscopic model for open quantum 
systems [9 considers the total Hamilton operator H = 
Hs + Hji + Hj^ composed of a Hamiltonian Hs of the 
system of interest, a reservoir part Hr and an interaction 
operator Hj of the form 
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Here, the system couples bilinearly to the reservoir via 
the dimensionless operator Q and the last term in i^/, 
the so-called counter term, guarantees that only the dy- 
namical impact of the reservoir is relevant. The dynam- 
ics of the density operator of the entire system W{t) 
obeys the Liouville-von-Neumann equation dtW(t) = 
—i/h[H^W{t)] with an initial state W(0) which factor- 
izes in a system part p(0) and the thermal equilibrium of 
the reservoir exp(— /3i^i^)/Z/3 at inverse temperature j3 = 

with partition function = TrR{exp(— 
Hence, the modeling of Hr in terms of harmonic oscil- 
lators assumes Gaussian noise properties of the forces 
i — Xln ^n^n of the heat bath with zero mean {(,{s))^ = 
and correlation function 
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mm)p = K'{s)+iK"{s) 



Jo 



coth 



cos(co'5) — i sin(co's) 



where the effective impact of the heat bath is captured 

by 

^^P^ = J^ I ds[K'is)[Q,[Qi-s),pit)]] 

+iK"{s)[Q,{Q{-s),p{t)}]}. (4) 

Here, the time dependence of Q{s) is taken in the in- 
teraction picture with respect to Hs and {• , •} denotes 
the usual anticommutator. The contribution ~ 7(0) 
in Eq. (|3| is related to the counter term in ([T]) with 
7(0) = fdw [JH/ttu;] 0. 

For the explicit representation of the operator master 
equation (|3| it is convenient to use the eigenstates of the 
bare system, i.e., Hs\n) = En\n). In the absence of ex- 
ternal time-dependent fields, the conventional procedure 
is to apply an additional RWA, where off-resonant, fast- 
oscillating terms are neglected. Thus, diagonal and off- 
diagonal elements of the density matrix pmn = {m\p\n) 
decouple to finally yield in the Schrodinger picture 

dp, 



(2) 



where the limit to a quasi-continuum of bath modes 
is performed by introducing a spectral density J{(jj) = 

f Xln m 1j ~ with dimension of energy. 

The relevant operator is the reduced density operator 
of the system, p{t) = TrR{W(t)}, which is obtained by 
tracing out the reservoir degrees of freedom. While this 
procedure does not in general lead to a simple equation 
of motion for progress is made in the weak-coupling 
limit upon employing the Born-Markov approximation. 
This yields the Redfield equation 



dp{t) 
dt 



,p{t) -jC[p]. 
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For the populations pnn{t) transitions between 
eigenstates are captured by rates Wnm = 
\{m\Q\n)\'^D{hjnm)/^'^, where Unm = {En - Em)/h 
denotes the bare transition frequency and 



D{huj) = 2Re 



ds K{s)e 



2hJ{uj)n(3{oj) (6) 



is the spectral noise density function with the single- 
particle Bose distribution n/3(cj) = l/[ex.-p{ujhp) — l]. The 
off-diagonal elements of the density matrix decay accord- 
ing to 

Tmn = ^=^^^±1^ + |M ((^IQI^) _ (^IQI^))' (7) 

with Tmm = ^k^m^km- The interaction between sys- 
tem and reservoir leads also to a shift in the bare tran- 
sition frequencies Comn = ^mn — ^^mn which is given in 
the Born-Markov limit by 
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where bath fluctuations appear in form of 

determined by 

K'^.^{oj) = / dsK'{s)sm{uJs) 
Jo 
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Above, the sine transformation 7sin and the Laplace 
transformation 7 of the classical damping kernel 7(5) = 
J duj [J{uj)/7ruj] cos{ujs) appear and Un = 27Tn/hP denote 
the Matsubara frequencies. The above expression for 
dujrnn capturcs the impact of reservoir fluctuations for 
both finite temperatures (ac Stark shift) and at zero tem- 
perature (Lamb shift). Since typical solid-state reser- 
voirs carry a broad range of spectral modes covering low 
as well as high frequencies, the system is embedded in 
an environment with a continuum of transition frequen- 
cies in contrast to the situation of a single-mode environ- 
ment such as a cavity. Accordingly, there are always low- 
frequency excitations present even at low temperatures 
and thus a clear separation between Stark and Lamb shift 
is not possible. Hence, in the sequel we simply use the 
notion Lamb shift for all reservoir-induced shifts of the 
transition frequencies. 

Let us specify the range of validity of the weak-coupling 
master equation ([5|. It applies as long as (i) a typical 
time for relaxation of the system towards thermal equilib- 
rium exceeds by far the retardation time of the reservoir 
and (ii) the level broadening induced by the heat bath is 
small compared to typical energy gaps of the bare sys- 
tem. Given a relaxation time of order I/7 (with a typical 
coupling strength 7 to the heat bath), retardation times 
of the bath of order 1/uJc (typical cut-off frequency cjc), 
and hp (thermal time scale) as well as a typical system 
frequency cjq, these conditions are fulfilled if 



jhP <C 1 and 7 <C cjo <C 



(11) 



This allows for very low temperatures compared with sys- 
tem energies such that the quantum optical limit j/uq <^ 
1 practically includes the zero-temperature limit. 



B. Two-level system 

To illustrate the above results let us consider a two- 
level system, namely. 



Hs = -Ez (7z + -Ex (7x , 



(12) 



where cr^ denote Pauli matrices and Ei are corresponding 
energies. In the experimental implementation to be dis- 
cussed below, the eigenstates of are associated with 
charge states which interact via Q = cf^ with a surround- 
ing reservoir so that [19^ 



Hi 



(13) 



Within the eigenstate representation with basis 
{|^), |e)} the system Hamiltonian (IT2]) reads 



(14) 



Hs = ^{\e){e\-\9){g\), 



where hwo = ^/E'^ -\- E"^. The operator Q takes the form 
Q = sm{vK\e){g\ + \g){e\) + cos(r?)(|e){e| - \g){g\) (15) 



with the mixing angle determined by tan(?7) = E^/Ez. 
This system is assumed to interact with a reservoir of 
Drude type 



J{uo) = ha- 



(16) 



where uOc ^ denotes a high cut-off frequency and a is 
a dimensionless coupling constant related to the typical 
dissipation strength discussed above via 7 aujQ. 

For the Lamb shift to be relevant, it must exceed the 
level broadening induced by the relaxation and excita- 
tion processes which is of order fiFge [see Eq. ([5|]. In the 
high temperature regime cjq^/^ <C 1, one has Fee ^ a/hf3 
and Eq. (|8| yields a Lamb shift that is negligible, namely, 
\SuJeg\ ~ (co'o/co'c) Fee- The rcservoir-iuduced Lamb shift 
is significant only at low temperatures uoh^ > 1, where 
Fee ~ c^ojo, while an explicit evaluation of Eq. ([8| to- 
gether with Eqs. ^ and (10) provides an effective tran- 
sition frequency uJe 



■ ujo — 5(jJeg with 



^ = ^Re Wiujo/J^i) - ^^{uJ,/u^)} . (17) 

Corrections to the above expression are of order ujq/uJc 
and t/j denotes the digamma function. For very low tem- 
peratures this result reduces in leading order to 
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which exceeds the level broadening by a characteristic 
logarithmic dependence on the cut-off frequency, also 
known from the atomic Lamb shift ll. For somewhat 
higher temperatures the Lamb shift ( |l7| ) decreases com- 
pared to the above asymptotic value. 



III. LAMB SHIFT OF A DAMPED HARMONIC 
OSCILLATOR 

In the case of a harmonic system, the reduced quantum 
dynamics can be calculated exactly for arbitrary spectral 
densities, temperatures, coupling strengths, and even in- 
cluding correlated initial states [9]. For factorizing initial 
states, an exact master equation with time-dependent co- 
efficients was first derived by Haake and Reibold |20l [21] . 
In the case of reservoirs with sufficiently large cut-off fre- 
quencies and in the regime jh^ < 1, these coefficients be- 
come time-independent after a transient period of time 
short compared to the relaxation dynamics. The cor- 
responding master equation allows to extract analytical 
expressions for the Lamb shift that can be used to spec- 
ify the range of validity of the weak-coupling results of 
Eqs. (Isl) and (17). In particular at sufficiently low tem- 



peratures the harmonic oscillator reduces to a two-level 
system [9 . Moreover, it provides insight in the non- 
perturbative regime of strong coupling which in turn will 
help us to identify the optimal parameter range to detect 
the Lamb shift in a broadband solid-state environment. 
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Hence, we study a system Hamiltonian Hs = p^/2M + 
Muj^q^ /2 of an oscillator with mass M which inter- 
acts via Q = Muoq/H with a reservoir of Drude type 
J{uo) = M'^uouo'l/ {uo'^ + ujD [note that for continuous de- 
grees of freedom it is convenient to use a spectral density 
which differs from that for discrete systems (16) by a 
factor Mujo/h with 7 = aujo]. One then obtains after a 
transient period for the reduced dynamics [20, 21^ 



Dp 



Q, [q, P{t)] 



D. 



p,[q,p{t)] . (19) 



The diflFusion coefficients take for uic ^» wq [H] the form 
Dp = 7(p^) and Dg = Mujl{q'^) - {p^)/M with the exact 
variances in thermal equihbrium [3 EQ] 
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(20) 



Upon representing Eq. ( |19| ) in the usual eigenstate ba- 
sis {|n)} of Hs and invoking the RWA, one derives in 
the subspace spanned by {|0), |1)} a transition frequency 
cDio = (jOq — 5uJiQ with 



E 



J7(k 



W0+'^n + k"l7(k"l) ' 



(21) 



Apparently, the Lamb shift is directly associated with the 
deviation from the equipartition relation, i.e., MujQ{q^) ^ 
(j9^)/M, induced by quantum fluctuations in the reser- 
voir. In agreement with the findings of the previous sec- 
tion, it is thus expected to be strong at very low tem- 
peratures. For weak dissipation, i.e., 7/cJo <C 1, the 
sum in Eq. ( [21] ) reduces in leading order to the sum in 
Eq. ( 10 ) which yields an expression of the form of Eq. ( 17 ) 
with a j/cjOq. We also note that the relaxation rate 
Tee ^ 7[(p^)/(^^^o) + 1] reduces in the weak-coupling 
limit to the known result Fge ~ 7[coth(cJo^/5/2) + 1]. 

The Lamb shift can also be extracted from Eq. ([T9| 
without applying the RWA. Restricting ourselves again 
to the subspace spanned by the ground and the first ex- 
cited state of the oscillator, the corresponding coupled 
equations for the density matrix elements can be easily 
solved to provide 



Wo 



2Dg 



(22) 



In the low-temperature range, an expansion for weak 
coupling gives duj'lf^^ ^ Dq/H + Dl/{2h^ujQ) + 



Dl/(2M^h^ujl) so that the RWA Lamb shift in Eq. (21) 
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FIG. 1: (Colour online) Lamb shift according to Eq. (21) for 
a Drude spectral density as a function of the dimensionless 
inverse temperature uoohp. The cut-off frequencies increase 
from bottom to top with uoc/ooo = 10 (dotted line), uoc/c^o — 
40 (dashed line), and uJc/ojq = 100 (solid line). The coupling 
constant \s ^/ujq — 0.05. The shaded region and the dots 
denote the temperature range and the Lamb shifts accessible 
for the discussed experiment. 



is obtained in leading order with corrections of order 
[(7/6^0 )ln(co'c/co'o)]^ at very low temperatures. Based on 
the full expressions, one can also show that in the range 
^hp < 1 the ratio |(5cj"^^^|/Fee decreases with increas- 
ing dissipation strength 7. This verifies the assumption 
that the low damping and high cut-off frequency domain 
constitutes the optimal detection window. 

We illustrate the above findings for typical parameters 
in realistic circuits to be discussed below. Figure [l] dis- 
plays the Lamb shift in the relevant weak-coupling case, 
for which Eq. (21) and the full expression (22) coincide 



up to minor deviations at low temperature. Whereas the 
shift depends on the coupling strength, the cut-off fre- 
quency, and the temperature, possibly the most conve- 
nient way to unambiguously identify the Lamb shift is to 
detect its temperature dependence. This way, an upward 
shift of about 6% of the bare transition frequency drops 
substantially towards higher temperatures. The optimal 
range of parameters for observing the temperature de- 
pendence of the shift corresponds to weak friction and 
large cut-off frequencies. 



IV. EXPERIMENTAL PROPOSAL 

In this section, we discuss a possible way to measure 
the Lamb shift for a heat bath with a broadband spectral 
density in a circuit containing a Cooper pair sluice [11 . A 
brief description of the actual set-up is given in the next 
section followed by an analysis of the detection proposal. 
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A. Physical system and its description 

The device we consider here, the sluice, is shown in 
Fig. |2] It consists of a superconducting island separated 
by two SQUIDs acting effectively as single tunnel junc- 
tions with tunable Josephson energies Jl.r- The electro- 
static potential on the island can be tuned by varying a 
gate voltage, Vg. The experimental access to the parame- 
ters Jl^r and Vg allows for a full control of the system and 
makes it an excellent prototype for detecting fundamen- 
tal quantum processes. Note that a similar device has 
already been used experimentally to study the geometric 
Berry phase related to adiabatic ground-state pumping 
of Cooper pairs [14| [22] . The Hamiltonian of the sluice 
can be expressed as [HI [15] 




Hs = Ec{N - Ngf - Jl cos (f - ^) - Jr cos (f + ^) , 

(23) 

where (p = (fn -\- (fL is the superconducting phase differ- 
ence across the sluice, Ng = CgVg/2e is the normalized 
gate charge, Ec = ^e^/Cj: is the charging energy of the 
sluice, Cg is the gate capacitance, and Cy: the total ca- 
pacitance of the island. Due to [Nj^^ (pj.] = —i the Cooper 
pair number operator Nk of the k-th SQUID has a rep- 
resentation A/'/e = —id^f^ {k = L^R), while the number 
operator of excess Cooper pairs on the island N is related 
to the corresponding phase difference 6 = {(pR — (Pl)/'^ 
hy N = —ide. If the device is operated in the charging 
regime, i.e., Ec ^ max{ J/,, Jr}, and the gate charge is 
close to a half integer, only the two lowest-lying charge 
states are important and we can adopt the two-state ap- 
proximation. Below, states with and 1 excess Cooper 
pairs on the island are denoted by |0) and respec- 
tively. 

This system is embedded in an electromagnetic envi- 
ronment which induces an effective transition frequency 
oj = uq — 6uj with uq and Su denoting the bare frequency 
of the two- level system (TLS) corresponding to Eq. (23) 
and the Lamb shift, respectively. Due to the broad 
bandwidth of the reservoir, the TLS cannot be isolated 
from its surroundings. This is in stark contrast, e.g., to 
a TLS interacting with a single-mode resonator studied 
in Ref. [4J. Our key idea is to identify the impact 
of the environment on the level splitting by adding 
an additional noise source consisting of a resistor R 
associated with equilibrium voltage fluctuations 5V as 
shown in Fig. [2j The resistor is directly coupled to the 
sluice island through a capacitor Ce- The advantage 
is two-fold: this reservoir can be engineered accord- 
ing to a given spectral density and the temperature 
of the sluice can be tuned by varying the tempera- 
ture of the resistor inducing a minimal effect on the 
rest of the circuit including the current detector, e.g., 
a large Josephson junction in parallel with the sluice [22] . 



The interaction Hamiltonian for the circuit in Fig. |2] 



FIG. 2: Circuit diagram for the Cooper pair sluice, the arti- 
ficial environment (in the dashed box), and the coil used for 
exciting the system with external drive (on the right). 



reads [ig [16] 



Hr 



-N2e^5V 
Cs 



(24) 



with voltage fluctuations determined by the fluctuation 
spectrum 



Rhuj 



1 + {cj/ujcy 



coth 



2kBT 



1 



(25) 



with the cut-off frequency uOc = {RCe)~^- The effective 
spectral density J{uj) of this reservoir is obtained from 
the spectral noise density in Eq. ^ at low temperatures, 
namely, limT^o D{— Huj) = 2hJ{uj). We find 



lim D(—huj) 



4e 



Ce 
Cy 

WH^ttR /Ce 



lim Sv(^) 



1 + {uj/ujcy 



,(26) 



where Rq = h/e^ is the resistance quantum. Hence, in 
the present circuit the reservoir induces in the relevant 
low-temperature regime a Drude-type spectral density 
(16) with the dimensionless system-environment cou- 
pling constant a = SttR/Rq {Ce/Cy)'^' Still within the 
range of validity of this approximation, we can effectively 
change the intensity of the Lamb shift by tuning the local 
temperature. 



B. Detection of the Lamb shift 

The possibility to control the parameters Jl,r and Vg 
enables the utilization of the optimal detection window 
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for measuring the Lamb shift. We fix the gate voltage to 
the degeneracy point A^^ = 1/2 and keep both SQUIDs 
open {Jl,r 7^ 0) but with different energies Jl Jr [23] . 
In this case, the Hamiltonian of the system (23) reduces 
to 

Hs{^) = -JLCO^[^-0)-jRCOs(^^^e) , (27) 

where in the restricted Hilbert space spanned by the 
charge states |0) and |1), the eigenstates are given by 



\9) 
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1 

71 



(|o) 
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(28) 



with fi = arctan j^qpj^ tan In this basis the Hamil- 
tonian is Hs{^) = {hwo/2){\e){e\ — \g){g\) with hwo = 
V^L + ^1 + 2 Jl Ji? cos (p. 

The basic idea is to analyse the spectrum of the two- 
level system, i.e., to measure its transition frequency, 
by exciting it using a microwave field with tunable fre- 
quency. The corresponding Rabi frequency is assumed to 
be sufficiently high such that the relaxation and dephas- 
ing processes do not play any role (see also the discussion 
at the end of Sec.|V|. Consequently, the interaction with 
the environment manifests itself only in a Lamb shift Suj 
which can formally be incorporated in the system Hamil- 
tonian Hs{^) ^ H^^^) = h/2 {uo - Suj){\e){e\ - \g){g\). 
A static phase bias (po is supplied by a magnetic fiux 
through the outer superconducting loop and a microwave 
field coupled inductively to the loop generates a time- 
periodic phase drive of the form S(p{t) = A cos[(cjo + ^^)^] 
with the off-set Sft around the bare frequency uq and 
small amplitude A <^ cpo. 

Let us study the driven system analytically before car- 
rying out the full numerical calculation in the next sec- 
tion. To this end, the total Hamiltonian Htot{^) ^ 
H^^{(po) + SHs is expanded up to the first order in Sep 
as 



SHs = [Jl sin ( Y - ^) + Jr sin + ^) 



With the help of the eigenstate basis in Eq. (28) it can 
be expressed as 



SHs = [hgg{\e){e\ - \g){g\) + /ige|ff)(e| + h.c] 



6(p 



(30) 



where hgg = J_ sin /i cos ^ — J+ cos /i sin ^ and hge = 



i{J- cos //cos ^ 



J+ sin /i sin ^ ) . 



The procedure is now straightforward. Upon switch- 
ing to a rotating frame with respect to (1/2) (cjq + t • 
(|e)(e| — \g){g\) in combination with a rotating wave ap- 
proximation one finds for the probability to reach the 
excited state at time t, provided the system is initially at 



t = prepared in the ground state, the known result 



"gel 



X sm 



t 

8h 



(31) 



Hence, in the off-resonance region + (^c^l ^ 

A\hge\/{4:h) the system is essentially trapped in the 
ground state and exhibits low amplitude, i.e., high- 
frequency oscillations, while in resonance SQ -\- Su ^ 
one has complete transitions with the Rabi frequency 

A\h,e\/{8h). 

The passage from the off-resonance to the resonance re- 
gion can be observed by monitoring the current through 
the circuit determined by the current operator / = 
2ed^pQHs/h since the ground and excited states support 
currents fiowing in opposite directions [25 . If the system 
is initialized in the ground state, a finite current fiows 
through the circuit. As soon as we enter the resonant re- 
gion, however, transitions are induced such that ground 
and excited state currents interfere destructively leading 
ideally to a vanishing net current in the circuit. The 
position of this current dip provides the total transition 
frequency of the system uq — Su. Its response to the 
temperature of the engineered environment gives access 
to the Lamb shift Suj. 



V. NUMERICAL RESULTS 

Here, we provide the full numerical simulation of 
the detection scheme for the Lamb shift introduced in 
Sec. 



IV 



The average net current through the sluice is 

presented in Fig. [3] as a function of the frequency off- 
set SQ. For the parameters employed here, we expect 
a current drop AI = (/max) — (/min) ^ 1 nA resulting 
in a detectable signal. By increasing the temperature 
from 0.05 K to 0.4 K, the resonance dip is shifted from 
0.063 X ujQ [solid curve in Fig. [3]; a)] to 0.033 X ujo [dashed 
curve in Fig. |3ja)] which is sufficiently large for being 
experimentally observable taking into account the typi- 
cal current noise at the detector [M]. The shape of the 
dip depends on both the current operator and the am- 
plitude of the induced Rabi oscillations, which can be 
varied, for example by tuning (po [see, for example, Eq. 
(31 )]. In Fig. [sjb), we show current dips for different val- 
ues of cpo. For decreasing cpo^ the current drop becomes 
sharper but the reference current decreases as well yield- 
ing a smaller difference AI. It turns out that for the 
optimal value (po = 7r/2, the dip is both sharp and deep. 
The experimental observation of the expected tempera- 
ture dependence of the frequency shift, given in Fig. [l] 
would provide conclusive evidence of the Lamb shift. 

Let us discuss the role of decoherence induced by the 
environment in our scheme. According to the general 




FIG. 3: (Colour online) (a) Average net current in a microwave- driven Cooper pair sluice as a function of the frequency shift 
6Q. The minimum current is found at SQ/uuq = 0.063 for T = 50 mK (solid line) and at 0.033 for T = 400 mK (dashed line). 
Currents are given in units of the critical current Ic = 3.14 nA yielding (/max) = 0.93 nA and (/min) = —0.06 nA. Data are 
averaged over a total time tfin = 7.6 ns with the system initialized to the ground state and the microwave drive turned on 
instantaneously at the beginning of the averaging period, (b) Same as in panel (a) but for various values of the phase bias 
^0 = 7r/2 (solid curve) and cpo = 0.5 (dashed curve) at T = 50 mK. The inset shows the normalized shape of the drop for the 
respective values of (po. Other parameters are Ec/ks = 1 K, A^^ = 1/2, Jl = 0.1 Ec, Jr = 0.05 /^c, and A = 0.2. 



theory presented in Sec. II A its impact is negligible as 
long as the frequency for Rabi oscillations is sufficiently 
higher t han th e decay rate of the coherences Tee ^ oo^Oi 
[see Sec. 
i.e., /^\ha 



II B 



in the relevant low temperature domain, 
, t'^ei ^ ^hw^a. To be consistent with this con- 
straint we have assumed a measurement time for the net 
current of t^^ = • 1000 ~ 7.6 ns which for typical 
device parameters still guarantees a unitary time evolu- 
tion of the sluice and in particular several tens of Rabi 
oscillations. 

Both the excitation of the system through the pertur- 
bation SHs and the measurement of the induced current 
must be done within tn^. In a realistic experiment, the 
main challenge in this situation is the synchronization be- 
tween the two pulses within tfin- If the pulses present an 
off-set, the average is less accurate since it is performed 
over a samller number of Rabi oscillations. However, se- 
quential measurements with different offsets should pro- 
duce a reliable current average. An alternative approach 
is to drive the system continuously, in the case of which 
there are no issues with timing the pulses and the mea- 
surement can be carried out slower. Floquet theory can 
be a good starting point to study this option in more 
detail [H]. 



VI. CONCLUSIONS 

We have presented a careful theoretical treatment of 
the Lamb shift including both the results for weak- 



coupling master equations and for a model of an exactly 
solvable harmonic degree of freedom. Contrarily to the 
situation of a cavity, where only a single-mode environ- 
ment is considered, we found it important to analyse the 
Lamb shift for thermal broadband reservoirs typical for 
mesoscopic solid state devices. Furthermore, a strategy 
to measure the Lamb shift in a realistic experimental set- 
up was discussed in detail supported by numerical simu- 
lations. Our results show a shift of about 6% of the bare 
transition frequency which is considerably larger than the 
so far measured Lamb shift of up to 1.4% of the qubit 
transition frequency in a superconducting waveguide res- 
onator [4 . In summary, the experimental realization of 
our proposal is well motivated. 
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